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In the MSSM with complex parameters, loop corrections to — * X^jW^ lead to 
a CP violating asymmetry A CP = (T(xt XjW+) - T(xT X°jW')) / (T{xf 
XjVF" 1 ") + — * XjW - ))- We calculate this asymmetry at full one-loop level. We 
perform a detailed numerical analysis for X\ ~* XiW^ and X2 ~* XiW^ analyzing 
the dependence on the parameters and phases involved. Asymmetries of several 
percent are obtained. We also discuss the feasability of measuring these asymmtries 
at LHC. 



It is well known that supersymmetric models contain new sources of CP violation if the 
parameters are complex. In the Minimal Supersymmetric Standard Model (MSSM), the 
U(l) and SU(2) gaugino mass parameters M 1 and M 2 , respectively, the higgsino mass 
parameter /i, as well as the trilinear couplings Af (corresponding to a fermion /) may 
be complex. Usually, M 2 is made real by redefining the fields. Non-vanishing phases of 
Mi and /i cause CP-violating effects already at tree-level in the chargino and neutralino 
production and decay [1, 2, 3]. In case the trilinear couplings of the third generation 
(A t , At,, A T ) are complex not only the stop, sbottom, and stau sectors [4] are strongly 
affected but also the Higgs sector [5, 6]. The three neutral Higgs bosons are no more CP 
eigenstates. 

Although new phases in addition to the CKM in the Standard Model (SM) are desirable 
to explain baryogenesis, there are severe constraints on the phase of \i from the exper- 
imental limits on the electric dipole moments (EDMs) of the electron, neutron and Hg. 
For example, in the constraint MSSM |0 M | has to be small [7, 8] for a SUSY particle 
spectrum of the order of a few TeV. 



In this note, we study CP violation in the decays xf(k 2 ) — > Xj(ki) + W + (—p) and 
Xi (— k 2 ) &i) + W~{p) in the MSSM with complex parameters by calculating 

the CP-violating asymmetry 

A _ r (+) fe + -> x°w + ) - r ( _)(xr -> x° 3 w-) 
CP r (+) (x+ -> %w+) + r ( _)(xr - x° 3 w-) 1 } 

at full one-loop order. The asymmetry is zero if CP is conserved and also vanishes at 
tree-level in case of CP violation. In Fig. 1 we show the graphs which contribute to this 
asymmetry at one-loop level. Of course, they give a contribution to Acp only if they have 
an absorptive part, i.e. some decay channels of xt must be open in addition to that into 

This asymmetry was already calculated in [9] considering only the third generation quarks 
and squarks in the vertex graphs. We have improved this calculation in several points. 
First, we performed a full one- loop calculation. In particular, we also calculated the 
contributions from self-energies of the charginos. It turns out, that these are important. 
(The self-energies of the neutralinos do not contribute due to their Majorana nature and 
the W ± - H ± transition vanishes for on-shell lU-bosons.) In addition, we take the Yukawa 
couplings running, which also gives a sizeable effect. Moreover, we take into account that 
the neutral Higgs bosons (h°, H°, A ) mix if the SUSY parameters mentioned are complex. 
In our case, this influence is, however, very small. As a loop-level quantity the asymmetry 
Acp depends on the phases of all complex parameters involved. One, however, expects 
that the dependence on the phases of Mi and A t)b is strongest (taking \i real). There is 
even a strong correlation between them. Therefore, a measurement of this asymmetry 
represents not only a test of CP violation in chargino decay, but can also be used for the 
determination of the phases of M 1 and A tjb . 
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Figure 1: All one- loop graphs of the decay xt ~ > Xj + W + , which contribute to the 
CP asymmetry Aqp defined in eq. (1), /' (/') denotes the isospin doublet partner of the 
fermion / (sfermion / ), e.g. f = t, f = b, <p° = (Hf, # 2 °, H°, G°), and 0+ = (H+, G+). 

The widths r (±) can be written as r (±) oc \Mlfi\ 2 + 2 Re[Mlf^M^ p ). 

The Feynman amplitudes for the tree- and one-loop level, jM^l and -M^,, are given by 

Mitel = 1^(^(0^ + L P L )u- t (k 2 )e;(-p), 

M£i = z%o(A; 1 )[(7 Ai Af + ) + mf +) )P R + (YAf +) + k^Uf +) )P L ]u^(k 2 )e;(-p) , 
MU P = ft^C-kJK^Afo + A^ (2) 

Since | A^t^el 2 — l-^^trel | 2 ? an d assuming, that the one-loop contribution is small compared 
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to the tree-level one, the CP-violating asymmetry Aqp takes the form 

ACP = \M^? ' (3) 

with the squared tree-level amplitude 

l-Mtreel 2 = P(\0 R \ 2 + \0 L \ 2 ) - 12m~+m,o Re[0 R *0 L ] , (4) 

and the one-loop contributions 

Re[M^M^ p ] = pRe[A R +) R * + Af +) L *] - Qm n m n Re[Af +) L * + Af +) R *\ 

+ 2^K° MK)° R * + n f + )° L i + m xt MK)° L * + K)° R *v > (5) 

Re[MUjMU p ] = pRe[A R _ ) R + Af_ } L ] - 6m~+m^o Re[Af_)0 L + Af^O R ] 

+^r( m x° M^f-)0 R + nf_ } o L ] + m-+ Re[nf_ } o L + nf_ } o fl ]) , (6) 

with the kinematic factor A = \(m 2 . + ,m 2 . ,my V ), X(x,y,z) = (x — y — z) 2 — 4yz, and 

,2 _ ™2 

-xj • -\° 

The chargino-neutralino-W coupling parameters L ' R , defined by the Lagrangian 



P = dt + 3 ( m l+ + m lo - m 2 w} 



c-wmt = w ^y(°fi p * + °% p L)xt + wfxti"{o R *p R + 0)*P L )x] , (7) 



are 



R = gZ* 2 U n + ^Z* 3 U i2 and 0% = gZ j2 V* - ^Z j4 V* , (8) 

where U, V, and Z are the matrices diagonalizing the chargino and neutralino system (see 
eqs. (A. 16) and (A. 17)). A and IT are form factors which are given in the Appendix A. We 
only give the form factors for x + and not for x~ , so that A, II always stands for A( + ), II( + ). 
The form factors A(_) and II(_), belonging to the x~ decay, can be easily obtained by 
conjugating all couplings. 



In Appendix A we present all formulas for the vertex contributions with the ttb and bbt 
loops and the chargino self-energy contribution with the tb loop, see graphs SF 1 F 2 , FSiS 2 , 
and SF of Fig. 1. The complete analytical formulas will be given in [10]. All individual 
one-loop graphs were numerically checked using the packages FeynArts, FormCalc, and 
LoopTools [11], and FF [12]. We included the CP-violating mixing of the neutral Higgs 
bosons by writing our own FeynArts model file. For the numerical program we used 
FeynHiggs [13]. 
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1 Numerical results 



We present numerical results for the decay rate asymmetries Aqp according to eq. (1) 
xf — > X^jW ± , for i = 1,2 and j = 1. A discussion of the other channels will be given 
in [10]. For the SM input parameters we take mz = 91.1875 GeV, my/ = 80.45 GeV, 
cos 6 W = m w /m z , tt{m z ) = 1/127.9, the on-shell parameters m t = 178 GeV, and m T = 
1.777 GeV. For the bottom mass, our input is the MS value mb(rrib) = 4.2 GeV. For 
the values of the Yukawa couplings of the third generation quarks (h t , hb), we take the 
running ones at the scale of the decaying particle mass. In principle, the parameters Af, 
the U(l) gaugino mass parameter M\ of the neutralino sector, and fi can be complex. 
We assume that \Mx\ = M 2 /2. In general, there are 15 independent sfermion breaking 
mass parameters. We take Mq as input and assume the MSUGRA inspired ratios rriq 



M, 



Q 



raj 



Mr 

M,2,3 



3:2:1 with m. 



Mr 



M T - 



M f 



Mb . 

^1,2,3 



In order to reduce the number of input parameters further, we use 



M T - 



M^ 3 , and 



A t = Ab = A T =: A. In all figures we take M40 = 300 GeV, tan/? = 10, and 0^ = 7r/10. 
For the decay 

Xi -> X?W /± , the total one-loop asymmetry Acp is shown in Fig. 2a and 
the tree-level branching ratio (BR) in Fig. 2b, for M 2 = 500 GeV, \A\ = 400 GeV, 
4>a = — 7r/4, 4>m 1 = 37r/4, and three values of Mq as a function of \fi\. \Acp\ increases 
for increasing values of because the tree-level decay width of Xi XiW ± goes to 
zero, as Xi becomes almost a pure bino which does not couple to W . Therefore, for 
Ia 4 ! ~ 550 GeV the branching ratio drops below 1%. The higher the value of Mq the 
heavier becomes the stop mass. Hence Aqp goes down but the branching ratio in (b) 
increases. 
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Figure 2: For M 2 = 500 GeV, \A\ = 400 GeV, (j) A = -vr/4, <j> Ml = 3vr/4, and three values 
of Mq, (a) the asymmetry Acp and (b) the tree- level branching ratio BR are given as 
functions of |/i|. 



Fig. 3 shows the dependence of A CP on (fi A for M 2 = 500 GeV, = 600 GeV, \A\ = 
400 GeV, and various 4>m 1 - A C p has its maximum at \4>a\ ~ tt/2 and is 



400 GeV, M Q 
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larger at large negative values of the phase <j>Mi- 
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Figure 3: The dependence of A CP on 0^, and various values of 4>Mii with M 2 = 500 GeV, 
|/x| = 600 GeV, |A| = 400 GeV, and Mq = 400 GeV. 

Now we discuss the asymmetry A CP for xt ~ * X\W ± . Fig. 4a shows the dependence of 
the asymmetry Aq P on the gaugino mass parameter M 2 for various values of \A\, = n, 
Mq = 300 GeV, <p A = -vr/4, and \fi\ = 200 GeV. For M 2 > 200 GeV, the lighter chargino 
and the two lighter neutralinos have dominating higgsino components and the heavier 
chargino is mostly gaugino-like (> 90%). The bigger \A\, the bigger is the mixing in the 
squark sector and hence Aq P . Around M 2 ~ 450 GeV the xt becomes massive enough 
so that the channels into bt 2 and £&i, 2 open. For M 2 > 250 GeV, the third generation 
(s) quark contributions clearly dominate the asymmetry, the self-energy contribution being 
bigger than the vertex contribution. For M 2 < 680 GeV, the vertex and the self-energy 
contributions for the third generation (s) quarks have opposite signs and cancel each other 
to a high degree. Nevertheless, they remain the dominant contributions in a large part of 
Fig. 4b. 

Various pseudothresholds are visible in Fig. 5a, where the squark mass parameter Mq is 
varied. The parameter set M 2 = 450 GeV, 4>\h = ^ an d I A* I = 200 GeV gives the masses 
= 468.55 GeV and m^o = 185.66 GeV. The strong dependence of A CP on the phase 
is clearly visible. Fig. 5b illustrates the dependence on the phase (f>A for Mq = {230, 300, 
400} GeV. That A CP does not factorize into a (pA dependent and a (pA independent part 
can be seen from the fact that the three curves do not meet in a single point. The other 
phases M and <pu x distort the factorization. 

The \n\ dependence of the decay xf ~^ Xi^ ± is shown in Fig. 6 for the same parameter 
set as used in Fig. 2. The total one-loop asymmetry Aqp is shown in Fig. 6a and the 
tree-level branching ratio (BR) in Fig. 6b. In the region ~ 400 GeV to 600 GeV the 
character of the xt an d Xi changes, for xt from gaugino to higgsino and for Xi from 
higgsino to mainly bino. Therefore, one has a strong dependence in Aq P and BR there. 
The dependence on Mq is analogous to that in Fig. 2. For > 600 GeV, the mass 
of xt ~ M an d Xi ~ M 2 /2 = 250 GeV. Therefore, the decay width of xt ~^ Xi^ ± 
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Figure 4: The dependence of v4 CP on M 2 for 
= 200 GeV, 

(a) the total asymmetry Aqp for various values of \A\. 

(b) the different contributions to the asymmetry for \A\ = 400 GeV: the vertex contri- 
bution with the third generation (s)quarks in the loop in black (1); the chargino self- 
energy contribution with the third generation (s)quarks in the loop in red (2); vertex and 
self-energy corrections with all other (s)fermions in the loop in blue (3); all remaining 
corrections in green (4). 
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Figure 5: For M 2 = 450 GeV, <p Ml =n,\A 

(a) the total asymmetry depending on Mq for various values of 4>a- 

(b) the total asymmetry depending on (pA for various values of Mq. 



increases with and Acp goes to zero. The hump in Fig. 6b at \y\ ~ 600 GeV for 
Mq = 500 GeV is due to the opening of the t 2 b channel. 

It is known that the electric dipole moments (EDM) of the electron, the neutron and mer- 
cury strongly depend on the phase of \i for a light SUSY spectrum [14]. The experimental 
constraints for the EDMs of the electron [15], the neutron [16], and mercury [17] can 
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Figure 6: For M 2 = 500 GeV, \A\ = 400 GeV, A = -7r/4, (j> Ml = 3vr/4, and three values 
of Mq, (a) the asymmetry Aqp and (b) the tree- level branching ratio BR are given as 
functions of 

be fulfilled by heavy sfermions of the first generations [18] or if cancellations of different 
contributions occur [8]. We checked for all plots all three EDMs and found always (small) 
values of <^ that fulfill all EDM constraints. 

Finally, we want to comment on the measurability of this asymmetry At LHC charginos 
are mainly produced in the cascade decays of gluinos and squarks so that the production 
rate strongly depends on their masses. If the gluino and squark masses are about the same, 
the gluino production cross section is far the dominant one. With rrtg ~ rriq = 750 GeV, 
we expect roughly 2.4 x 10 5 events containing xf (one has the same amount of xt an d Xi 
in the case where they originate from gluinos or from a gluon-gluon process), assuming 
a luminosity of 10 5 pb _1 and a branching ratio of a gluino decaying into a xt °f 40%. 
Taking into account the branching ratio for xf — >• Xi^ 1 ^ one can measure Aqp for this 
decay with a statistical significance of ~ 2 (confidence level of 95%). For measuring Aqp 
for xf ~^ Xi^ 1 ^ assuming 5 x 10 4 events containing a xt or X21 one g e ^s a similar 
statistical significance. 



2 Conclusions 

We have calculated the CP- violating asymmetry between the partial decay rates T(xf — > 
X^W + ) and T(x~ — > x°jW~) due to phases in the MSSM. It is a pure loop effect. We 
have calculated this asymmetry at full one-loop order. We have given numerical results 
for Xi — * Xi an d xf — * Xi^ ± - The respective asymmetries are of the order of several 
percent, depending on the values of parameters and phases involved. In order to have 
reasonable branching ratios for the decays the Xi rnust not be very bino like. We also 
discussed the feasability of measuring such an asymmetry at LHC. It might be possible 
to measured it with a confidence level of 95%. 
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A The (s)top/(s)bottom contributions to Acp 



The relevant parts of the Lagrangian are 

£ = -9=w;lrfP L t - J^W+t-fPLb (9) 
-\-ic mn W^ b m d^t n + ic mri W^ t n d^b m 

(jmi^R ^mi-fi) Xi b m Xi (jmi^L + k mi Pjij t b~ m 
+b (fniPR + KnPl) X~t n + l~ (tP L + fe£Pfl) bt* n 

+i(ai J P R + b l n] P L ) xf n + I] (a%P L + b%P R ) tt* n 
+b (a b mj P R + b b mj P L ) x% + X°j (a%P L + b^Pn) b b* m , (10) 
with the coupling parameters to the W boson, 

C mn = ~~^J^mL^nL J (H) 



(12) 



(13) 



the chargino, 

4 = -gVaKl + ghtV l2 K\ , k l ni = g h b u* 2 K* L , 
L = -gUiA + gKu i2 vt R , k h mi = g h t v*nl; L , 

and the neutralino, 

<• = gKlflj - gh t K* R z j4 , b% = - g h t ni* L z* 4 + g ni* R f Rj , 
<• = gKtJh - ghK; R z j3 , b b m3 = - g h b n b n : L z*, + g ii b * R f% 

with the gaugino components of the neutralino 

f% = V2e q t^6 w Z* l , p Lo = -V2({e q -II L )t^e w Z ]1 + II L Z ]2 ) , (14) 
and the Yukawa couplings 

\J2m w sin/3 1 \/2in w cosp V / 

The charge and the isospin of the quark q are given by e q and 7f L , g is the SU (2) coupling 
parameter. 

The mixing matrices are defined as 

UjcMtfVfr = 5 jk m^ , (16) 
Z- a M%Zip = S jk m^o for the basis {6, ^ 3 ,^,^} , (17) 

jaV vl Ja/3'^k/3 ~ °jk m2 q k 



nUM 2 )Lni = 5 jk ml , (1 
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with the mass matrices 
M* + = 

M*° = 

(M 2 y = 



m w cos/3 pi J 

/Mi -m z s w cp m z s w Sf3 \ 

M 2 m z c w Cf3 -m z c w s j3 



-m z s w c f3 m z c w Cf3 
V m z s w sp -m z c w sp 





-fj, 



-n 





/ 



m~ L a q m q 
a q m g m\ R 



where the abbreviations sw = sin^, cw = cos 9w, sp = sin/5, cp = cos/3, and 
m\ L = M| + m|cos2/3(/| L - e(? sin 2 ^)+- 2 



m 



\ R = M~ R + e q m| cos 2(3 sin 2 6 W + 



in 



A q - (j,* (tan (3) = \a q \e 



o i( t>q 



(19) 
(20) 
(21) 

(22) 

(23) 
(24) 



are introduced for a more convenient notation. M~ R = M^iM 2 ^) is the soft breaking mass 
parameter for the right stops (sbottoms). 

The form factors giving the major contribution for most of the parameter regions studied, 
can be split into the contributions (ttb), (bbt), and {lb): 



\ L,R a L,R , a L,R , a L,R 

t-jL^R _ T-rL,R | t-tL,R 

n (+) - n (ttb) + L \bbt) ' 



with 



IV (ttb) 



Iy (ttb) 



LL {ttb) 



11 (ttb) 



iy {bbt) 



•^-^= E a^m^o(A;**(Co + Ci)m 6 + Z£(C + Ci + C 2 )m*+) 
+ ^m t (^C m 6 + C(C + C 2 )m~ + ) , 



-^j= 2 E 46(^0 - 2C 00 + + dmlj + C 2 m 2 , 



) 



+ ^njlniCim^mt + a t nj k t * i C 2 m b m^+ , 



E <.(^C 2 m b + Z£(C 2 + Ci 2 + C 22 )m-+) 

3 2 _ _ 



g 2 Z-^ u mj fh ni u mn , ~ y 00 j 



(25) 
(26) 



(27) 

(28) 

(29) 
(30) 
(31) 



m,n=l 
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3 2 



tfbbi) ~ 0Z2 £ ab mjlni C mnC 00 , (32) 



\bbt) g 2 

m,n=l 

3 2 



n fe = ^E c+ n (<.e,(C 2 + C12 + C 22 K+ + fc^-A^Cd + C 11 + C 12 )m^o (33) 



m,n=l 

- a* -4i(C + Ci + C 2 )m 6 ) , 



2 



ngsij = A E ^(a^-e^ + C u + C 12 )m^o + &^(^ + C 12 + C 22 )m-+ (34) 

°" m,n=l 

- 6^£(C + d + C 2 )m 6 ) , 
A & = ~^ -2 £ MV.)^ + ^(H^ M ) (35) 

+ (So + S 1 )m-+(4 (3 _ i) ^>-+ + ^(3- i )^>x+_ i) ) > 

x « X (3-i) 



+ (S + J Bi)m- + (/J l(3 _ j) ^>- + + k^k^m- 



(3-i)' 



where 



Cj* b) = C x (m|o, m 2 ^, m 2 + , m? ,m 2 ,m^) , (37) 
B? b) = B (m 2 w ,m 2 t ,m 2 b ) , (38) 



C ( x bt) = C x {m\^m 2 w ,m\+, m 2 , m\ ,m~ ) , (39) 
,(&) 



i?r = B x (m^,ml,mi) . (40) 
The B- and C-functions are given in the notation of [19]. 
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